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Abstract. Quantum walks on the line with a single particle possess a classical 
analog. Involving more walkers opens up the possibility to study collective quantum 
effects, such as many particle correlations. In this context, entangled initial states and 
indistinguishability of the particles play a role. We consider directional correlations 
between two particles performing a quantum walk on a line. For non-interacting 
particles we find analytic asymptotic expressions and give the limits of directional 
correlations. We show that introducing ^-interaction between the particles, one can 
exceed the limits for non-interacting particles. 
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1. Introduction 



Quantum walks were introduced []J as a generalization of a classical random walk to 
a unitary evolution of a quantum particle. The time evolution can be either discrete 
or continuous jij. The connection between discrete-time and continuous-time quantum 
walks has been established for a walk on a line j^l, 0] and, more recently, for walks on 
arbitrary graphs Q. It has been shown that both continuous js| and discrete time 
quantum walks can be regarded as a universal computational primitive. Continuous 
time quantum walks have been extensively studied in the context of coherent energy 
transfer in networks 10]. Both continuous and discrete time quantum walks have found 



a promising application in designing quantum algorithms 
algorithms based on quantum walks have been proposed 12 



Indeed, a number of 



21 , for a review see 22 



Various properties of quantum walks have been analyzed, in particular their 



asymptotic behavior 



23 



-25 



sec 



and the effect of the initial conditions 



26-28 



for a review 



29]. Due to the wave-nature of quantum walks a number of counter- intuitive 



ohenomena has been observed, including infinite hitting times 30J, |3l| and localization 



321-138]. The properties of random walks on infinite regular lattices are closely related to 
the dimensionality of the lattice. It is well known that a classical random walk returns 
to the origin with certainty in dimension one and two while in higher dimensions the 
probability of return (Polya number) is strictly less than unity 39]. For the discrete- 



time quantum walk the recurrence properties are determined not only by the dimension 
but also by the initial state and the coin operator leading to rich behavior 40h45| . The 



46 



closely related property of persistence has been studied in 

The extensive theoretical studies have stimulated the search for experimental 
implementations of quantum walks. Various schemes based on ion traps 471 ]. optical 



lattices [48j, |49j, cavity quantum electrodynamics [50|, optical cavities [5l| or Bose 



Einstein condensate 52j have been proposed. Recently, discrete time quantum walk on 



the line has been realized in a variety of physical systems including cold atoms 
trapped ions [54j, |55[ and photons [56j, [57 . 



53 



Most of the studies to date considered quantum walks with a single particle. A 
natural extension of the field of quantum walks is to involve more particles. This 
unlocks the additional features offered by quantum mechanics such as entanglement 
and indistinguishability which are not available in classical random walks. Quantum 
walk on a line with two entangled particles has been introduced in 58] and the meeting 
problem in this model has been analyzed 59]. A physical implementation of this model 
based on linear optics has been proposed |6Q_J. Quantum walks with two particles have 
been applied to the graph isomorphism problem 6jJ, [62j . Entanglement generation in a 

Recently, 



63 



special two-particle quantum walk on a line has been investigated in 
the first successful experiment with two particles on a line has been reported 64 
A framework for multi-particle quantum walks on rather arbitrary graphs has been 



proposed in 65]. The study of quantum walks with more particles on the line is 



motivated by the fact that the single-particle walk in this case can be considered as 
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a classical interference phenomenon [66J. We note that walks on higher- dimensional 
lattices cannot be considered classical in this sense, since the resources needed to 
simulate the quantum walk scale exponentially. 

In this paper, we investigate the non-classical effects in the two-particle discrete- 
time quantum walk on the line by asking the question: How is the directional correlation 
affected by the quantum nature of the particles? In particular, we analyze the probability 
P s of finding both particles on the same (negative or positive) half-line. We derive 
analytical expressions for the asymptotic value of this probability in dependence on 
the initial coin state. Classically, a symmetric random walk has a fixed value of the 
probability P s equal to 1/2. We first consider two quantum particles on a line starting 
the walk in a separable state. We determine the limits for the directional correlations and 
show that, for any value within these limits, one can design a corresponding separable 
initial state. Next, we prove that the bounds cannot be exceeded by considering 
entanglement in the initial state. On the other hand, introducing quantum walks with 
5-interactions, we show that the directional correlations can be increased above the limit 
for non-interacting particles. 

Our paper is organized as follows: we briefly review the quantum walk on a line 
with one and two non-interacting particles in Section [2] and introduce the probability 
to be on the same side of the lattice P s . In Section [3] we analyze the probability P s 
for separable initial states. Entangled initial states are considered in Section HI In 
Section [5] we study the influence of the indistinguishability on the probability P s . In 
Section |6] we introduce the concept of ^-interacting quantum walks to break the limits 
of non-interacting quantum walks. We summarize our results in Section [7J 

2. Quantum walk on a line with one and two particles 

Let us first briefly review the quantum walk of a single particle on a line (see e.g. Ref. 



67( for a more detailed introduction). The Hilbert space of the particle is given by a 
tensor product 

U = U P ®Uc 

of the position space 

Up = £ 2 (Z) = Span{|m)| m G Z} 

and the two-dimensional coin space 

He = Span {\L),\R)} . 

We consider a particle starting the quantum walk from the origin, i.e. the initial state 
has the form 

|V(o)> = |o)®hk7>, 

where \ipc) denotes the initial state of the coin. After t steps of the quantum walk the 
state of the particle is given by 

|^(t)> = (Mm,t)\m)\L)+Mm,t)\m)\R)) = £7*^(0)) , (1) 
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where the unitary propagator U has the form 

U = S(I®C). (2) 

The coin operator C flips the state of the coin before the particle is displaced. In 
principle, C can be an arbitrary unitary operation on the coin space He- We choose the 
most studied case of the Hadamard coin, denoted by Ch, which is defined by its action 
on the basis states 

Ch\L) = ^=(\L) + C H \R) = ±=(\L) - \R}). 

After the coin flip the step operator S displaces the particle from its current position 
according to its coin state 

S\m)\L) — ► \m-l)\L), S\m)\R) — > \m + l)\R). 

The coefficients ipL(fy(m,t) in (JTJ represent the probability amplitudes of finding the 
particle at position m after t steps of the quantum walk with the coin state \L(R)). The 
probability distribution generated by the quantum walk is given by 

p(m,t) = \(m\(L\m)\ 2 + \H(R\m)\ 2 = \M™,t)\ 2 + \M™,t)\ 2 . 
The extension of the formalism described above to two distinguishable particles has 



been given in 58] . One should consider the bipartite Hilbert state as a tensor product 

Hu = H X ®H 2 

of the single particle Hilbert spaces. We consider non-interacting particles, i.e. their 
time evolution is independent. Hence, the propagator of the two-particle quantum walk 
can be written in a factorized form 

U 12 = U 1 ®U 2 , (3) 

where U\ (U 2 ) is the propagator of the first (second) particle given by Eq. (j2j). Note 
that this factorized time evolution cannot increase entanglement between the particles. 
In this paper we consider particles starting from the same lattice point (the origin). 
Hence, the initial state of the two-particle quantum walk has the shape 

|*(0)) = |0,0)<g>|tt o ), 

where \^c) is the initial coin state of the two particles. 

Let us first consider the case when the initial coin state is separable, i.e. 

|*c> = ® ■ (4) 

Since entanglement is not induced in the process of time evolution, the two-particle state 
remains factorized and the joint probability distribution p(m, n, t) of finding the first 
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particle at the mth and the second at the nth site at time t is reduced to the product 
of single particle distributions 

p(m,n,t) =pi(m,t) -p 2 (n,t) . (5) 

Here, Pi(m,t) is the probability distribution of a single-particle quantum walk given 
that the initial coin state was Hence, the two-particle quantum walk with initially 
separable coin state is fully determined by the single-particle quantum walk. 

We turn to the situation when the initial coin state \^c) does not factorize. In such 
a case, the joint probability distribution p(m, n, t) cannot be written in a product form 
([5]). Nevertheless, we can map the two-particle walk on a line to a quantum walk of a 
single particle on a square lattice. Indeed, we can write the two-particle propagator (jHJ) 
in the following form 

U l2 = S X2 {I Pl2 ®{C H ®C H )), (6) 

where Jp 12 is the identity on the joint position space and the joint step operator S\ 2 
is given by the tensor product of the single particle step operators S^. The relation 
(J6]) implies that we can consider the two-particle propagator U\ 2 as a propagator 
of single-particle walk on a plane with the coin given by the tensor product of two 
Hadamard operators. Hence, the two quantum walks in consideration are equivalent. 
This correspondence allows us to treat the joint probability distribution of the two- 
particle walk with the tools developed for the single-particle quantum walks. 

Finally, let us briefly comment on a quantum walk with indistinguishable particles. 
It is natural to use the second quantization formalism. We denote the bosonic creation 
operators by at and the fermionic creation operators by bt ~, e.g. at ■•. creates one 
bosonic particle at position m with the internal state \i), i = L, R. The dynamics of the 
quantum walk with indistinguishable particles is defined on a one-particle level, i.e. a 
single step is given by the following transformation of the creation operators 



o} n — >-Lfa in +of +1 J, a\ m -+J-(a] in -aJ +1 J 

(m,L) fK V [m— l,L) (m+l,R) J ' (m,R) /K \ (m—l,L) (m+l,R) J 



for bosonic particles, similarly for fermions. The difference is that the bosonic operators 
fulfill the commutation relations 



0, 



"K«)'"(nj) 



$mnO~ij i (7) 



while the fermionic operators satisfy the anti-commutation relations 

b( m ,i), b(nj) }=0, tyn,j) } = ^rnn^ij . (8) 

Since the dynamics is defined on a single-particle level, one can describe the state of 
the two indistinguishable particles after t steps of the quantum walk in terms of the 



single-particle probability amplitudes (see Ref. [59| for a more detailed discussion). 

In the present paper we focus on the directional correlations between the two 
particles. We quantify this property by the probability P s that both particles are found 
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after t steps of the quantum walk on the same side of the line. For distinguishable 
particles it is given by 

oo t t 

m=~t n=—t ra=l n=l 

For indistinguishable particles p(m,n,t) = p(n,m,t), i.e. these two probabilities 
correspond to the same physical event. Hence, the sums in ([9]) have to be restricted 
over an ordered pair (m, n) with m > n, i.e. 

n=— t \m=n / n=l \m=n / 

In particular, we will be interested in the asymptotic limits of the probability P s in its 
dependence on the initial coin state of the two particles. We consider both separable 
and entangled coin states, as well as indistinguishability of the particles, in the following 
Sections. 



3. Separable initial states 

Let us now specify the probability P s for two distinguishable particles which start the 
quantum walk with a separable coin state (j4j. As discussed in the previous Section the 
joint probability distribution p(m,n,t) factorizes (JSJ). Therefore, the probability to be 
on the same side of the lattice P s simplifies into 

P s (t) = Pf(t) • P 2 (t) + P+(t) ■ P+(t) , (11) 

Here we have denoted by P i (i) the probability that the particle which have started the 
quantum walk with the coin state is on the positive or negative half-axis after t 
steps, i.e. 

i 
m=—t m=l 

In Figure [1] we plot the course of the probability P s (t) with the number of steps 
t. To unravel the dependence on the initial coin state \^c) we consider three cases 
- (i) = \L) ® \R) (black dots), (ii) \^ c ) = \L) ® \L) (open circles), and 

(Hi) \^ c ) = \{\L) +i\R)) ® ^(\L) +i\R)) (open diamonds). We find that after 
some initial oscillations the probability P s quickly approach steady values which are 
determined by the initial coin state. 

Let us now determine the asymptotic value of the probability P s in dependence of 
the initial coin state. Consider a general separable coin state of the form 

\^ c ) = (ai\L) + h\R)) ® (a 2 \L) + b 2 \R)) . 



The asymptotic probability distribution for a single particle is given by 24 



p(x, t, a u bi ) = ^f^ + ^ + ^-W . (12) 
Wl-2S(l-#) 
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Figure 1. The probability P s that two distinguishable particles performing a quantum 
walk on a line end on the same side as a function of time. Both particles start the 
quantum walk from the origin. As the initial coin state \^c) we choose one of the three 
factorized states - (i) \L) for the first particle and \R) for the second particle (black 
dots), («) \L) for both particles (open circles), and (m) \ips) = -7= (\L) +i\R)) for 
both particles (open diamonds). We find that for the initial coin state (i) the particles 
are more likely to be on the opposite side, since P s < 1/2. Indeed, due to the choice of 
the coin state \LR) the probability distribution of the first particle is biased to the left 
while the probability distribution of the second particle is biased to the right. On the 
other hand, for the initial state \LL) both probability distributions are biased to the 
left. Hence, the particles are more likely to be found on the same side. Finally, for the 
initial state (Hi) which results in the symmetric single-particle probability distribution 
the particles are equally likely to be on the same or the opposite side of the line. The 
asymptotic values of P s for all three initial states are in agreement with the analytic 
estimation of Eq. (fT5l) . 



The probability that the particle is on the negative or positive half-axis is obtained by 
integrating the probability density over the corresponding interval 



P l ~(a i ,b i ) = J p(x,t,ai,bi)dx = ^(2 + ((a i + b i )ar i + (ai~bi)bi)) , 
t_ 

V2 
t 

P^(a h bi)= J p(x,t, ai,bi)dx = ^(2 - ((a* + bi)a~ + (a^ - &*)&»)) . (13) 


Note that within the approximation of Eq. (TP2]) the resulting integrals are time- 
independent, i.e. we immediately obtain the asymptotic values of the probabilities 
Pj. This is due to the fact that the asymptotic probability density depends only on 
the ratio x/t. 
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Inserting the results of (1131) into the Eq. (Hip we find that the probability P s is given 

by 

p(se P ) = 1(4+ (( ai + h)al + (ax - h)h) ((a 2 + 62)02" + (a 2 - 62)62") ) . 

(14) 

In particular, for the initial states (i — iii) considered in Figured! we find the asymptotic 
values 

Pj L/? ) = P s (l, 0, 0, 1) = I , Pj LL) = P s (l, 0, 1, 0) = I , 

p ™ ^ Ps( 71'71'7!'71 ) = ^ (15) 

These results are in perfect agreement with the numerical simulations presented in 
Figure [H 

Let us now analyze the probability P^ sep > in more detail. First, we recast the formula 
( fl^l) in a simpler form by a change of the basis of the coin space. Consider the basis 
formed by the eigenstates of the Hadamard coin 

c H \x ± ) = ±\x ± ), 

which have the following expression in the standard basis 



W-) = ^\L)±^\R). (16) 
We decompose the single-particle coin state in the Hadamard basis 

m = K\x + ) + K\ x -). 

From the expression ( Tl6l) we find the transformation between the coefficients in the 
standard and the Hadamard basis 



-hf + - hj , h = - hj - hj . 



2 2 ' ' ' 2 2 

With the help of these relations we find that the formula (1141) for the probability P^ sepS> 
simplifies in the Hadamard basis into 

Pi^) = i(2 + (2|/ i +| 2 -l)(2|^| 2 -l)) . (17) 

Here we have used the normalization of the single-particle coin state \ipi), i.e. the 
condition 

\ht\ 2 + \K\* = l. (18) 

We display the probability to be on the same side pj 56 ^ in its dependence on 
the parameters hf in Figure |2j We find that pj sep - ) reaches the maximum value 3/4 
provided that both hf equals zero or unity, i.e. when both particles start the walk in 
the same eigenstate of the Hadamard coin. Indeed, starting the single-particle walk in 
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Figure 2. The probability to be on the same side P^ Bep ' in its dependence on the 
coefficients of the initial coin states. The parameters hf are given by the overlap of 
the coin state \tpi) with the eigenstate |x + ) of the Hadamard coin. We find that the 
probability pj sep ^ reaches the maximum value 3/4 when both h\ 2 equals either zero 
or one. The minimum value 1/4 is obtained if one the hf is zero while the other one 
is unity. 



the eigenstate \x + ) {\x )) leads to a probability distribution which is maximally biased 
towards left (right). We illustrate this feature in Figure [3) Note that this effect has been 



identified numerically in [26]. Hence, when both particles start the walk in the same 
eigenstate of the Hadamard coin, their probability distributions are maximally biased 
towards the same direction and, consequently, the particles are the most likely to be 
on the same side. On the other hand, if the particles start the walk in the different 
eigenstates (e.g. the first one in \x + ) and the second one in |x~), which corresponds to 
— 1 and h^ = 0), the probability distributions are maximally biased in the opposite 
directions. In such a case, the particles are the most likely to be on the opposite side of 
the lattice and P^ sep ' reaches the minimum 1/4. 

4. Entangled initial states 

Let us now analyze the probability that the particles will be on the same side of the 
lattice P s for the initial coin states \^c) which are not factorized. We follow two 
approaches. First, we analyze the particular case of maximally entangled Bell states. 
We express the two-particle state in terms of single-particle amplitudes. In this way, 
we decompose the joint probability distribution into single-particle distributions plus an 
interference term. We then use the results of the previous section to find the asymptotic 
value of the probability P s . By this approach we emphasize the role of the interference 
of probability amplitudes. Second, we employ the equivalence between the two-particle 
walk on a line and single-particle walk on a square lattice discussed in Section |2j This 
correspondence allows us to use the tools developed for the quantum walks with a 
single particle, namely the weak limit theorems 23|, to find the asymptotic probability 
density for the two-particle walk on a line. We leave the details of the calculation for 
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Figure 3. Single-particle probability distribution for the initial coin state \i/ic) — \x + )- 
We find only one peak on the left side of the lattice, the peak on the right side has 
disappeared. Consequently, the resulting probability distribution is maximally biased 
towards left. Choosing the initial coin state as \ipc) = |x~) will &ip the plot around 
the origin and the resulting probability distribution will be maximally biased to the 
right. 



(\LL) ± \RR)) 



(19) 



the Appendix A With the explicit form of the probability density we finally derive the 
asymptotic value of the probability P s for an arbitrary two-particle coin state. 
We start by examining the particular case of maximally entangled Bell states 

±(\LR)±\RL)) , 1^) = ^ 

Obviously, the joint probability distribution p(m, n, t) is no longer a product of the 
single-particle probability distributions. However, we can still express it in terms of the 
single-particle probability amplitudes. Let us denote by ip^ (to, t) the amplitude of the 
particle being after t steps at the position m with the coin state \i), i = L, R, provided 
that the initial coin state was \L). Similarly, let ip\ (to, t) be the amplitude for the 
initial coin state \R). With this notation we express the joint probability distributions 
generated by quantum walk of two particles with initially entangled coins by 

p^ (to, n, t) = J \^ L) ( m ' ] ( n > *) ± ^i R) ( m ' *)^f } fa t) | 2 > 

i,j=L,R 

\4 L \m,t)^ L) (n,t) ±4 R \m,t)i;f\n,t)\ 2 , 

~ i,j=L,R 

where the superscript indicates the initial coin state. We now make use of the fact that 

(T R\ 

the amplitudes ip\ ' (to, t) are real valued. Indeed, both the Hadamard coin and the 
initial states have only real entries. Hence, the amplitudes cannot attain any imaginary 
part during the time evolution. Therefore, we can replace the absolute values by simple 
brackets and expand the joint probability distributions in the form 

p(^ ± )( 77i) n ,t) = - 



p^ \m, n, t) 



- £ 

i,3=L,l 



t)p {R \n,t) + p( R \m,t)p( L \n,t) ) ± 
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± ip(m, t)<p(n, t) , 
p (0±) (m,n,t) = - (p iL) (m,t)p {L \n,t) + p (R) (m,t)p {R) (n,t) ) ± 

±<p(m,t)<p(n,t). (20) 
Here, we have used the notation 

to shorten the formulas. When we insert the expressions (1201) into the definition (Q of 
the probability P s we find that the later one can be written in the form 

p(^)(t) = P( LR \t) ± 1(f), P^\t) = P { s LL \t) ± lit) . 

The interference term I(t) is given by 

I(*)=(^(t))%(W)) 2 , 

where we have denoted 

o t 
<p~(t) = ^ (p(m,t), <^ + {t) = ^2^(m,t) . 

m=—t m=l 

Let us now turn to the asymptotic values of P s in dependence on the choice of the 
Bell state. The limits of pj LR ^ and pj LL ^ are given in (fT5T) . We obtain the asymptotic 
value of the interference term I(t) from the numerical simulation, which indicates 

J(f^+oo) = i. 

Finally, for the limiting values of the probability P s we find 

pW>+) = I ; pW>") = I ? p(4> + ) = ^ pW") = I _ (21) 

We display the dependence of P s on the number of steps and the choice of the initial 
coin state in Figure HI We find that the probability P s quickly approach the steady 
values, similarly as for the factorized coin states which we have shown in Figure [H For 
(open circles) and \<fi~) (black dots) the particles are asymptotically equally likely 
to be on the same or on the opposite side. For the Bell state (stars) the particles are 
more likely to be on the same side of the line. Finally, for the singlet state \if)~) (open 
diamonds) the particles are more likely to be on the opposite side. The asymptotic 
values of the probabilities P s are in agreement with the results of ( |2TI) . 

After we have analyzed the particular case of the Bell states we turn to a general 
initial coin state. As in the previous Section, we make use of the asymptotic probability 
density p(x\, X2, t) and replace the sums in ((9]) by integrals. We derive the explicit form 



of the asymptotic probability density in the Appendix A Performing the integrations 
we arrive at the following expression 



P s = ^2 + \h {++) \ 2 + \h 



|2 17 |2 It |2 
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Figure 4. The probability that two distinguishable particles performing a quantum 
walk on a line with initially entangled coins end on the same side as a function of 
time. Both particles start the quantum walk from the origin. As the initial coin state 
l^c) we choose one of the Bell states (TIT))) . For the Bell states \ip + ) (open circles) and 
\4>~) (black dots) the particles are equally likely to be found on the same or on the 
opposite side of the line in the long time limit. For finite times they are more likely 
to be on the same side for \tp + ) and more likely on the opposite side for \<fi~). For the 
other two Bell states (open diamonds) and \(f> + ) (stars) the differences remain in 
the asymptotic limit. The particles are more likely to be on the opposite side for the 
singlet state \tp~) and more likely to be on the same side for \<fi + ). The asymptotic 
values of the probability P s agree with the findings of (f2Tj) . 



for the probability to be on the same side. Here we have denoted by h( a p) the coefficients 
of the decomposition of the initial coin state \^c) into the basis formed by the tensor 
product of the eigenvectors \x k ) of the Hadamard coin C#, i.e. 

\*c) = E h i^)\x a )\x P ). (22) 
Finally, using the normalization condition for the initial coin state \^c) 

\h {++) \ 2 + \h h+) \ 2 + |/i (+ _ } | 2 + |/i ( __)| 2 = 1, 
we can simplify the expression for the probability P^ ent ^ into the form 

n (ent) = ^(l + 2(|/.( ++ )| 2 + IV-)| 2 )) • (23) 

The dependence of the probability pj ent - ) on the initial coin state is illustrated in 
Figure [5j We find that the probability to be on the same side for entangled initial coin 
states P^ ent ^ satisfies exactly the same bounds as the probability p^ sep ^ derived in the 
previous Section for separable initial coin states. The maximum value of 3/4 is reached 

when |/i( ++ )| 2 + \h/ )| 2 = 1. In such a case, the initial coin state \^c) is an eigenstate 

of the two-particle coin Ch <8> C# corresponding to the eigenvalue +1. On the other 
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Figure 5. The probability to be on the same side of the lattice Pj in its 
dependence on the choice of the initial coin state. We find that _pj en *- > is bounded 
in the same way as pj sep ^ displayed in Figure [21 The maximum is obtained for states 
satisfying the condition |/i(++)| 2 + \hi )| 2 = 1, while the minimum is reached when 

= h—) = °- 



hand, the minimum value 1/4 of the probability P s is attained when both 7ir++) and 

/i( ) vanishes. This corresponds to \^c) being the eigenstate of the coin Ch®Ch with 

the eigenvalue —1. 

Finally, we note that for separable coin states the formula ( 1231) reduces to Eq. (TlTl) 
which we have derived in the previous Section. Indeed, for separable states we have the 
relation 

which together with the normalization ( fl8|) implies 
P s = I^l + 2 |/ i +| 2 |/ i +| 2 + 2|^| 2 |/i 2 | 2 ) 

= l(l + 2|^| 2 |4| 2 + 2(l-|^| 2 )(l-|^| 2 )) 
= i(2 + (2|^| 2 -l)(2|4| 2 -l))=Pi^). 

5. Indistinguishable particles 

Let us now briefly discuss the probability to be on the same side P s for indistinguishable 
particles. We show that for a particular choice of the initial state of the two bosons or 
fermions the problem reduces to the case of distinguishable particles with maximally 
entangled coins. 

As the initial state of the quantum walk we choose 

|*(0)>=|1(0,L)W, 

i.e. both particles are initially at the origin with the opposite coin states. Recalling 
the amplitudes ip\ L ^ ) for the single particle performing the quantum walk with 
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the initial coin state \L) (\R)) we express the state of two bosons and fermions in the 
following form 

m,n i,j=L,R 

i* (F) w> = E E ^('"^^(Mii^UiH, (24) 

m,n i,j=L,R 

where \vac) denotes the vacuum state. Note that in (jUj) both summation indexes m 
and n run over all possible sites. Using the commutation (JJJ) and anti-commutation (JHJ) 
relations we can restrict the sums in ( 124")) over an ordered pair (m, n) with m > n. The 
resulting wave-function will be symmetric or antisymmetric. 

We turn to the joint probabilities p(m, n, t) that after t steps we detect a particle 
at site m and simultaneously a particle at site n, with m > n. First, for m/nwe find 



P 



(m,n,t) = I <l (m>4) l (nJ) |*( s .- p ) (t)) I 



= E \4 L \rn,t)^f\n,t)±4 R \m,t)^ L \n,t) , 

i,j=L,R 

where the + sign on the right hand side corresponds to the bosonic (B) , and the — 
sign to the fermionic (F). Comparing these expressions with the results for Bell states 
( )20l) we identify the relation 

p {B \m,n,t) = 2p^ +) (m,n,t) , p (F \m, n, t) = 2p { ^\m, n, t) . (25) 

For m = n we obtain for bosons 

p( B \m,m,t)= |(2 (miL) |^ B )(t))| 2 + |(2 {m , fl) |vI/( B )(t))| 2 + 



+ 1(1 



{m,L)i-{m,R) 



¥ B \t))\ 



= 2 



2 



.(*)/ 



+ 



+ 



(m, t)^f > (m, t) + # (m, t ) 4* ) (m, t) 
and for fermions 

p^(m,m,t) = \(l {m , L) l {m , R) \¥ F \t))\ 2 = ^ L \m,t)^\m,t)-^\m,t)4 R \m,t) 



We note that relations similar to ( |25l) hold as well for m = n. Indeed, we find the 
following for bosons 

p( B \m,m,t) = - |^( L ^(m, t)tpj R \m, t) + x/j^ R \m, t)ipj L \m, t) 

i,j=L,R 

= pW + \m,m,t) , (26) 

and for fermions 

p( F \m,m,t) = - \^ L \m,t)^pj R \m,t) — ^ R \m,t)^pj L \m,t) 



i,j=L,R 



(27) 
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Finally, we derive the probability P s that the bosons (or fermions) are on the same 
side of the line. As we have already discussed, for indistinguishable particles we have 
used the formula ( flOi) where the summation is restricted to an ordered pair (m, n) with 
m > n. However, using the results of ( 125|) . ( 126|) and ([27)1 we can replace p( B,F \m,n,t) 
by p^ \m,n,t) in (flOl) and extend the summation over all pairs of m and n. Hence, 
we find that 

pw(t) = p« & \t), pW(t) = p«-\t). 

In summary, the results for bosons (resp. fermions) are the same as for distinguishable 
particle which have started the quantum walk with entangled coin state (resp. 

This is a direct consequence, of course, of the required symmetry properties of 
two-particle boson and fermion states. We note that also the fact that the particles 
have started the walk from the same lattice point is important. However, when 
the two indistinguishable particles start the walk spatially separated their evolution 



59 



Indeed. 



differs from that of distinguishable particles with entangled coin states 
indistinguishability starts to play a role when the wave-functions begin to overlap, 
whereas entanglement is a non-local property. 

6. Quantum walks with 5-interactions 

We have seen in the preceding sections that entanglement in two-particle non-interacting 
quantum walks cannot break the limit of probabilities we found for separable particles. 
A natural question arises: What happens if we consider interacting particles? This 
motivates us to introduce the concept of two-particle quantum walks with ^-interaction. 
To do that, we change the factorized time evolution operator defined in (|3]). In the 
original time evolution the coin was the same factorized coin in all lattice point pairs 
(m, n), in the 5-interaction quantum walk we change the coin to a non-factorized one 
C<5, when the particles are at the same lattice point m = n. 

Considering the above, we define the unitary time evolution operator for quantum 
walks with ^-interacting particles on a line as 

U s = S 12 (P S ® (C H ® C H )) + S 12 (P S ® C s ) , 

where P$ is the projector on the joint position state 

Ps = \m)\m)(m\(m\ , 

m 

and 

Ps = Ip 12 — Ps ■ 

As an example, we consider the entangling 5-interaction coin Cg of the following form 

/ 1 1 1 l\ 

1 -1-11 



-11-11 
V-i -i i iy 



(2f 
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Figure 6. Joint probability distribution (upper plot) and the probability to be on 
the same side of the lattice P s (lower plot) for two interacting particles performing a 
quantum walk on a line. The (^-interaction coin C$ is realized by a matrix (|28jl. As 
the initial coin state we have chosen one of the Bell states, namely \^c) — The 
resulting joint probability distribution is mostly concentrated on the diagonal, as can 
be seen from the upper plot. Consequently, the particles are very likely to be on the 
same side of the lattice. Indeed, the lower plot indicates that the asymptotic value of 
the probability P s exceeds 0.8. 



In Figure [6] we present the results of a numerical simulation of the corresponding 
quantum walk with ^-interaction. The initial coin state was chosen to be the Bell state 
From the upper plot we find that the joint probability distribution is concentrated 
on the diagonal, thus the particles are likely to be found on the same side. The lower 
plot clearly indicate that quantum walks on a line with 5-interactions can break the 
upper limit of P s = 3/4 which we have derived for non- interacting particles. 
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7. Conclusions 

We have analyzed the two-particle quantum walk on a line focusing on the directional 
correlations between the particles. The directional correlation of two non-interacting 
particles on the line is shown to be confined in an interval, independent of wether the 
initial state is entangled or not. The bounds of the interval are reached when the initial 
states coincide with the eigenstates of the coin operator. 

Introducing a ^-interaction one can exceed the limit we derived for non-interacting 
particles. The 5-interaction breaks the translational symmetry, thus new analytical 
tools are needed to investigate the properties of the introduced model. In the picture 
of the joint time evolution, this scheme could be considered as an inhomogeneous two- 
dimensional quantum walk, where the coin is changed on the diagonal line m = n. 

Acknowledgments 

The financial support by MSM 6840770039, MSMT LC 06002 and the Czech-Hungarian 
cooperation project (KONTAKT,CZ-11/2009) is gratefully acknowledged. SMB thanks 
the Royal Society and the Wolfson Foundation for their support. 

Appendix A. Asymptotic probability distribution for a quantum walk with 
two entangled particles 

In this appendix we derive the asymptotic probability density for a quantum walk on a 
line with two particles for an arbitrary initial coin state \^c)- We make use of the close 
relation between the two-particle walk on a line and a single-particle walk on a plane 



23 



discussed in Section [21 We then employ the weak limit theorem 

The time-evolution of the Hadamard walk on a plane is in the Fourier representation 
determined by the propagator 

U 12 (k u k 2 ) = U 1 (k 1 )®U 2 (k 2 ). 

Here, Uj(k) denotes the single-particle propagator of the Hadamard walk on a line, 
which is given by 

U J (k)=V(e- tk ,e lk )-C H . 

Since U\ 2 (ki, k 2 ) has a structure of a tensor product of two unitary matrices we write 
its eigenvalues in the form 

Atf(Jfei, k 2 ) = e iuJ ^ klM) = e ^(^)+^fe)) ; z, j = 1, 2 , (A.l) 

where e luJi ^ are the eigenvalues of the matrix Uj(k). Their phases Ui(k) are determined 

by 

(sin k \ 
— j=- J , u 2 (k) = 7r — Ui(k) . (A. 2) 
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Similarly, we write the corresponding eigenvectors of U\ 2 {ki, k 2 ) i n the form of a tensor 
product 

Vij(h, k 2 ) = Vi(k x ) ® Vj(k 2 ) 
of the eigenvectors of the matrices Uj(kj) 



vi(k) 



L= (e ik , V2e iuJl{k) -e ik ) 



\/ni(k) 



1 T 

v 2 (k) = - (-e ik , V2e' iu)l{k) + e lk ) . (A.3) 
Vn 2 (k) V / 

The normalization of the eigenvectors is given by 

ni(k) = 2 ^1 + cos 2 k — cos ky/l + cos 2 k^ , 

n 2 {k) = 2(1 + cos 2 k + cos kVl + cos 2 k^ . 
The weak limit theorem 23[ states that the cumulative distribution function equals 



F(Z u x 2 )=J2 [ 



JJ= 1 ^ Vw i j 1 ((-oo,xi)x(-oo,a;2) ) 

where we have denoted Xj = y. The probability measure fiij is determined by 



(A.4) 



Vij = \(v i j(k 1 ,k 2 ),^c)\ 



2 dki dk 2 
2vr 2tt 



The four-component vector corresponds to the initial state of the coin l^c 1 ). From 
the explicit form of the eigenvectors Vij(ki, k 2 ) given in (I A. 3 1) we find that the probability 
measures fi^ equal 

1 [ i + (~iy +1 ( y c 1 c(k 1 ) + s 1 s(k 1 ^ + 
(-iy +1 (c 2 c(k 2 ) + s 2 s(k 2 ) 



(A.5) 



(-l) l+J ( C 12 C(h)C(k 2 ) + 5 12< S(A; 1 )«S(A: 2 )+ 

dfci dk 2 
2tt 2tt ' 



+X 1 C(/c 1 )5(/c 2 ) + X 2 S{kx)C{k 2 
Here, we have used the notation 

cos k sin 



C(fc) 



VI + cos 2 fc V1 + cos 2 /c 

to shorten the formulas. The coefficients C, S and X entering the expressions flA.5j) 
can be determined from the initial state of the coin \^c)- 

To obtain the cumulative distribution function (1A.4|) we also have to find the 
integration domains. These are determined by the gradients of the phases Uij(ki,k 2 ) 
of the eigenvalues of the propagator Ui 2 (ki, k 2 ). From their explicit form given in (1A.1I) 
and (IA.2j) we find that the gradients are 

VuyikxM) = ((-l) m C(A;i), (-iy +1 C(k 2 )) . 
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Using the above derived results and the substitution 

riu\ coski ai d ^ 

C{ki) = = = q h dki 



we can simplify the cumulative distribution function into the form 



F(xi,x 2 ) 



X\ X2 

1 f dqi f dq2 



1 - C x q x - C 2 q 2 + C 12 q x q 2 



x/2 

With the help of the relation 



d 2 F 

p(x,y) 



dxdy 

between the cumulative distribution F(x, y) and the probability density p(x, y) we find 
that the later one is given by 



vr 2 (l - d) J I - 2#(1 - 4) J I ~ 2# 



Xx n X 2 _ X\X 2 

- Cxj - C 2 - + C 12 — 

Finally, we give the explicit form of the coefficients Cx, C 2 and Cx 2 - We find that they 
have a particularly simple form in the basis formed by the tensor product of eigenvectors 
of the Hadamard coin Ix^ 1 ), which have been given in f|T6|) . With the decomposition 
of the initial coin state in the Hadamard basis as given in (|22p we obtain the following 
expressions for the coefficients Ci j2 and C\ 2 . 

Cx =>/2(|/ H++) | 2 + |/i (+ _ ) | 2 -|/ l( _ +) | 2 -|fc ( „ ) | a ) , 
C 2 = v ^(|/ l(++) | 2 + |/, ( _ +) | 2 -|/ l(+ _ ) | 2 -|/, ( __ ) | 2 ) , 

c 12 = 2 (i^i 2 + 1 v-)l 2 - l^ (+ -)l 2 - 1 V+)l 2 ) • 
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